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NONPARAMETRIC  LITE  TEST  SAMPLING  PLANS 

1.  Introduction 

Sampling  plans  for  truncated  life  teste  based  on  several  parametric 
families  of  distribution  functions  have  been  proposed  by  Sobel  and 
Tischendorf  (1959),  Gupta  and  Groll  (1961)  and  Gupta  (1962*)-  Thus, 
we  have  the  problem  of  which  family  to  chooee;  i.e.,  exponential, 
normal,  lognormal,  gamma,  etc.  In  many  cases  it  is  assumed  that  certain 
nuisance  parameters  are  known.  Unless  considerable  prior  information  is 
available,  the  choice  will  necessarily  be  somewhat  arbitrary.  Ve  can 
circumvent  this  problem  by  considering  a  class  of  distributions  restricted 
only  by  certain  intuitive  considerations.  If  the  items  in  question  wear 
out  with  age,  then  it  is  reasonable  to  assume  that  the  failure  rate 
function  is  nondecreasing.  For  other  items  such  as  solid  state  electronic 
devices,  it  seems  reasonable  to  assume  that  the  failure  rate  function 
is  nonincreasing.  By  using  sharp  bounds  on  such  distributions,  we 
can  obtain  sampling  plans  which,  although  more  conservative  than  the 
parametric  plans,  offer  greater  protection  since  they  are  valid  for  a 
much  larger  family  of  distributions. 

It  has  been  suggested  by  Zelen  and  Dannemiller  (1959)  that  sampling 
plans  should  be  written  in  terms  of  percentiles.  Such  plans  were  presented 
by  Gupta  and  Groll  (1961)  and  Gupta  ( 1962a)  for  the  gamma,  normal  and  log¬ 
normal  distributions.  We  present  plans  based  on  percentiles  under  the 
increasing  (decreasing)  failure  rate  assumption  as  well  as  plans  based 
on  the  mean  life.  A  partial  discussion  of  these  plans  was  also  pre¬ 
sented  by  Gupta  (1962b). 


2.  Properties  of  Increasing  (Decreasing)  Failure  Rate  Distributions. 


Let  X  be  a  random  variable  with  distribution  F  such  that 
F(x)  =  0  for  x<0.  If  F  has  density  f  then  r(t)  = 
f(t)/[l  -  F(t)j  is  touvn  as  the  failure  rate  function.  Note  that 
r(t)  ■  -  ~  log  (1  -  F(t)j  when  a  density  exists.  For  this  reason, 
we  say  that  F  is  UR  (DFR)  for  increasing  (decreasing)  failure 
rate  if  log  (1  -  F(t)]  is  concave  where  finite  (convex  on  [0  ,  «>)). 

These  distributions  have  been  studied  by  Barlow,  Marshall,  and  Proschan 
(19b3)  and  Bariev  and  Marshall  (1964).  For  convenience  we  sunmarize 
same  of  the  fundamental  properties  of  these  distributions. 

An  UR  (DFR)  distribution  ray  have  at  most  one  point  of  dis¬ 
continuity;  namely,  at  the  end  of  its  right  (left)  hand  interval  of 
support.  The  continuous  part  of  an  UR  (EFR)  distribution  is 
absolutely  continuous.  Bums  of  independent  1FT*  random  variables  theiwelves  have 
UR  distributions.  This  is  not  true  for  EFR  random  variables.  If 
F  is  a  mixture  of  EFR  ’istributions,  then  it  follows  that  F  is  EFR. 

Many  of  the  properties  of  the  UR  (DFR)  family  f ollcw  from  comparison 
with  the  exponential  distribution.  This  is  not  surprising  since  it  is 

the  boundary  distribution  between  the  two  classes.  If  U,  <  U_  <  . . .  U 

1  2  n 

are  the  order  statistics  from  an  UR  ^R)  distribution  F  with  mean 

n  . 

n  then,  E[U.]£h.  while  E[U  )  £  u  Z  ±  . 

^  1  (^)  n  (£)  1  1  * 

To  obtain  sampling  plans  based  on  percentiles  we  use  the  fact 

that  [1  -  F ( t )  ] ^  ^  is  decreasing  ( increasing)  in  t  when  F  is 
UR  (DFR).  Let  p.  denote  the  r^  moment  of  F  .  To  obtain  sampling 

plan/  based  on  bounds  on  moments  we  use  the  inequalities  f  r  UR  dis¬ 
tributions.  Now  we  describe  the  relevant  inequalities. 
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(2.1)  1  -  F(t)  £  < 


j"  exp  [  -  t/Ty 

Lo 


-Vxi/r)  ,  — 1/r 


(r  £  1) 


,  t  >  n; 


l/r 


(2.2)  1  -  F(t)  £  < 


r 

L 


t  <;  ^/r 


(r  >  0) 


exp  [  -vt]  ,  t  >  m.; 


l/r 


where  w  uniquely  satisfies 


.  r  -  1  wx. 
lx  *  r  i  x  e  dx 


and 


\  -  nr/r(r  +  1) 


th 


If  F  is  EFR,  with  r  moment  ^  , 


then 


(2.3)  1  -  F(t)  £  < 


exp  [  -  t/>yr]  ,  t  £  r>^  r 


r  -r 

re  Hj. 

T(r  ♦  i)tr 


,  <■  rX^/r 


th 

From  these  bounds  we  can  obtain  the  following  bound  on  the  q 
quantile  £  In  terms  of  ^  assuming  F  is  IFR.  If  q  £  1  -  e 
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then 


[  -  log  (1  -  qHuj  S.  C,  £  1  "  log  ^  ^  i 


if  q  £  1  -  e 


1 


Frca  the  fact  that  X^r  is  decreasing  (increasing)  in  r  when 


The  density  and  the  failure  rate  are  hounded  at  the  origin  frc*  above 
in  teras  of  acaents  as  follows: 


f(0)  £ 

Ai+J-1 


1,  2, 


when  F  is  UR. 

Bounds  on  densities  and  failure  rates  vere  obtained  by  Barlow 
and  Marshall  [196U].  If  F  is  UR  and  r  £  1  then 


f(t*)  £  *t+)  £  ^ 

*V  -  t 


1/r 


°  ^  t  < 


0  ,  t  ^  ^ 
f  ^  r 

•(*->  ^  l  l/r 

1  *  ,  t  >  iC 


when  a  satisfies  4 


f*  r  -  1 
.  -  r  /  x  1 

r  J0 


-ax. 
e  dx 


If  F  is  ETR  vlth  density  f  ,  then 


(tej 


-1 


s  <V 


1/r 


f(t  )  ^  <  X/'  e 


l/r  -t/V1/r 


.  *J/r s»s(r* D^/r 


a! 

Xr(£-~)  e"(r4l)  ,  t  £  (r  ♦  l)X^/r 


Finally,  ve  mention  that  maximum  likelihood  estimates  of  IFR  ( DFR) 
distributions  have  been  obtained  by  Grenander  (1956)  and  Marshall  and 
Pros chan  (1964). 


A 
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3-  Nonparamctrlc  sampling  Plans 

Let  ©  denote  a  parameter  defined  by  a  life  distribution  F  . 

th 

For  example,  0  might  denote  the  mean  life  or  the  90  percentile. 

We  would  like  to  establish  that,  say  ©  ©^  where  ©^  is  some 

specified  value.  A  cannon  practice  in  life  testing  is  to  truncate 

the  experiment  at  a  preassigned  time,  say  t  ,  and  note  the  number  of 

failures  (assuming  that  a  failure  is  well-defined).  One  object  of  these 

experiments  is  to  set  a  confidence  (lower)  limit  on  the  mep.n  (quantile)  life  of  the 

units.  If  It  is  desired  to  establish  a  mean  (quantile  life,  with  a  given  proba- 

* 

bility  of  at  least  P  and  if  the  number  of  failux'es  in  the  fixed  time 
does  not  exceed  a  given  number  c  ,  then  the  experimenter  would  like  to 
know  the  minimum  sample  size  necessary  to  achieve  his  objective.  In 
following  these  sampling  plans  and  the  associated  decision  rule  the 

probability  of  rejecting  a  lot  having  a  parameter  ©  £  ©^  is  at  ltaot 

* 

equal  to  P  Alternatively,  the  consumer's  risk  in  adopting  the  sampling 

* 

plan  (n  ,  c  ,  t  ,  ©  )  does  not  exceed  1  -  P  whatever  ©  may  be. 

#  * 

Mathematically,  given  a  number  P  (0  P  £  1)  ,  a  time  t  and  a  value 

©^  of  ©  and  an  acceptance  number  c  ,  we  want  to  find  the  smallest 

positive  integer  n  so  that  if  the  observed  number  of  failures  in  time 

# 

t  does  not  exceed  c  we  can  assert  with  confidence  probability  P 
that  ©  ;>  ©  . 

It  should  be  pointed  out  that  these  life  tests  can  be  termin.ted 
prior  to  the  time  t  with  rejection  as  the  result.  In  fact,  the  ter¬ 
mination  of  the  experiment  car.  take  place  at  the  smaller  of  the  two 

s  t 

times  t  ,  t,  ,  \  where  t,  .  is  the  time  to  the  (c  ♦  1)  failure. 

’  ( c+1)  (c+I ) 

Also,  the  associated  decision  rule  leads  to  an  acceptance  only  at  the 
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end  of  tine  t  and  only  if 


*(0*1) 


>  t 


Hence, 


L(p)  •  P( Acceptance)  *  ^l“(c+l)  ^  ^ 

(5.1)  -  ;c  .  Tjfin  -  c):J  -  f(t;#]n'c‘1ar(t;8) 

t 

C 

Z/n\  i,_  Nn-i 

(j)  P  (1  -  P) 

1-0 

where  p  -  F(t;G)  . 

If  p  Is  b  decreasing  function  of  ©  which  is  true,  for  instance, 
if  the  density  f(t;0)  is  TP^Totally  positive  of  order  2)  in  t  and 
©  ,  then  F(t;0)  £  F(  t;0Q)  <— >  ©  £  ©Q  .  If  pQ  -  F(t;©0)  ,  then  it 
is  clear  from  the  theory  of  confidence  intervals  that  the  smallest 
to  satisfy 


c 


i-0 


is  such  that  in  adopting  the  above  sampling  plans  with  the  associated 

* 

decision  rule,  we  can  assert  with  probability  P  that  ©  £  ©Q  provided 
*(0+1)  >  *  • 

It  should  be  noted  that  the  above  discussion  assumes  a  knowledge 
of  F(t;0)  .  New,  we  consider  the  case  where  F(t;0)  is  not  explicitly 
assumed;  it  is  only  assumed  that  we  have  bounds  on  7(t;0)  . 

Suppose  that  p  -  F(t;0)  £  b(t;©)  for  t  £  0  where  b(t;©)  is 
a  known  function.  Suppose  further  that  b(t;0)  is  decreasing  in  ©  . 
Since  L(p)  is  a  decreasing  function  of  p  ,  we  have 
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(5-5) 


L(  P  ( t ;  © ) )  £  L(b(t;0) )  £  L(b(t;©0)) 


for  ©  £  ©^  .  Hence  the  sampling  plans  of  the  above  type  are  obtained 
by  choosing  the  smallest  positive  integer  n  satisfying 

(3-4)  L(b(t;©0))  £  1  -  P* 

where  c  ,  ©^  and  t  are  fixed 

From  (3-3)  and  (3-4)  it  follows  that  if  the  number  of  failures 

during  time  t  is  less  than  or  equal  to  c  ,  then  we  can  make  the 

* 

confidence  statement  that  ©  £  ©^  with  confidence  probability  P 

Note  that  if  b(t;©)  is  a  sharp  bound  on  F(t;9)  ,  then  it  is 
nondecreasing  in  t  since  F(t;0)  is  nondecreasing  in  t  .  If  © 
is  the  mean,  then  b(t;©)  =  b(t/©;l)  since  the  mean  is  a  scale  parameter 
for  positive  random  variables  and  b(t;©)  is  decreasing  in  ©  . 

Now  we  give  a  few  examples  where  the  bound  b(t;©)  is  known. 

Example  1.  (Markov's  inequality).  Let  ©  denote  the  mean.  Then 

0  ,  t  <  © 

F(t;0)  £  b(t;0)  -  \ 

L  1  -  ©/t  ,  t  ^  * 


is  true  for  nonnegative  random  variables.  Clearly,  b(t;©)  is  de¬ 
creasing  in  9  .  A  sampling  plan  based  on  Markov's  inequality  would 
afford  protection  over  the  class  of  all  distributions  on  the  positive 
axl6.  However,  the  bound  is  quite  wide. 
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Example  2.  (Unimodal  density) 


If  the  density  f  is  unimodal  with  unknown  mode  and  first  moment 
©  and  If  F(O-jQ)  -  0  ,  then  from  Barlcw  and  Marshall  (1965),  we  have 


2© 


F(t;©)  £  b(t;0)  «  <2-~  , 


t<4 


1  - 


2t 


Clearly,  b(t;G)  Is  decreasing  in  0  and  the  bound  is  a  slight 
Improvement  over  Example  1. 

Example  3.  (PF2  density). 

If  log  f ( t ;©)  is  concave  where  f(t;©)  >  0  ,  then  f  Is  called 
a  PFg  ( P<5lyn  Frequency  of  Order  2)  density.  Most  of  the  coon  only  used 
densities  nave  this  property.  If  f  has  mean  0  ,  then 


F(t;0)  ^  b( 


t;0)  =  1 

L  1  - 


»  *  £® 

sup  (1  -  e  bt)/(l 
n£>t 


-bm  v 

•  )  »  t  ^  , 


•/; 


xbe’bxdx/  (1 


-bm. 

e  )  • 


where  b  is  determined  by  © 

The  above  bound  Is  tabulated  in  a  paper  by  Barlow  and  Marshall  (I96J). 
A  sampling  plan  based  on  this  bound  would  afford  less  protection  than 
the  UR  sampling  plane  since  f  PF^  implies  that  F  is  EFR  but  not 


conversely. 
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for  the  UR  and  QFR  Distributions 


New  we  discuss  the  derivation  of  sampling  plans  for  the  UR  and 
DFR  distributions  for  the  following  two  cases. 


UR  Distributions 


Case  (  i  )  W  =  u. 

r 

In  this  case,  from  (2.2),  we  have 


(3-5) 


t » M-  )  =  ^ 


^  1  -  e-Vt 


*  <s  K 


t  >  u 


where  w  uniquely  satisfies 


(3-b) 


r-1  -wxJ 

=  r  /  x  e  dx  . 


Clearly  b(t;p  >  is  a  decreasing  function  of  .  Hence  the  required 
n  to  insure  >  p^  is  given  by  the  smallest  positive  integer  which 

satisfies 


(3-D 


L  J  ‘"''Ur  ,  U 


b(t;u  ;  J  •  S  1  ’  P 
r  - 


It  should  be  noted  that  the  solution  of  (3-7)  depends  on  the  fact  that 
r  0 

t  >  p,  .  ‘Of  specia^  importance  is  the  case  where  r  =  1  .  In  this 
case  (3.6)  reduces  to 


(3-0)  41 


w  =  1  -  e 


Hence,  using  (3-d)  to  solve  for  w  and  then  using  (3-7)  with 
b(t;p^)  =  b(t/p^;l)  =  1  -  e  Vt  ,  sampling  plans  to  insure  a  specified 
mean  life  have  been  c  mputed  and  are  given  in  Table  1.  This  table 
also  gives  the  actual  confidence  level  attained. 


Case 


(u)  -  -  c. 

In  this  case 


the 


th 

q 


percentile. 


(3-9) 


*  <  c. 


which  is  again  decreasing  in  £  and  clearly  depends  only  on  q  and 

q 

the  ratio  t/£  .  Hence,  to  insure  a  specified  4uantile  life  ^  , 
we  substitute  from  (3-9)  the  value  of  b(t;£°)  in  the  general  equation 
(3.4)  and  solve  for  the  smallest  n  .  Table  II  at  the  end  of  this 
paper  gives  minimum  sample  sizes  for  the  life  testing  situation  where 


DFR  Distributions 

Case  ( 1 )  W  »  p  . 

r 

In  this  case,  letting  >  ■  p, ,^T(r  ♦  1)  ,  ve  have 


1  - 


(3-10)  b(t;p  ) 


^  1  -  rr> 


.(et) 


-r 


t  <  r>4  r 

t  i  r)y/r 


and  one  can  construct  sampling  plans  for 

o1' 1 

and  t  £  r(>  )  .  For  the  speclSLl  case 


0 

the  cases  where  t  £  r  (>  ) 
r  »  1  ,  euid  t  <  p^  ,  it  is 
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easy  to  check  that  the  sampling  plans  baaed  on  the  bounds  in  ('.1J) 

reduce  to  the  sampling  plans  f  >r  the  negative  exponential  with 

0  0 

A^  =  p.-^  =  ^  These  tables  have  already  been  given  by  Jobel  and 
Tischendorf  ( 1 /;  >  )  for  t  <  >  =  u 

Case  ( i i )  w  =  ^  . 

In  this  case  the  bound  is 


VC 

1  -  U  -  q)  q  ,  t  ^  cq 

(3.11)  b  ( t ,  C  )  *  I 

q  1  0  .  t>?q  . 

To  use  tiie  jpomtrivdjlU.‘bound,  we  require  t  <(  £ ~  .  The  sampling  plane 

are  obtained  as  in  case  (.11)  of  IFR  distributions  discussed  above. 

Table  III  at  the  end  of  this  paper  gives  these  sampling  plans  for 
*  ■  t/Cq  i  1  • 

PC  Functi un_  and  the  Producer's  Risk 

The  UC  function  represents  the  probability  of  acceptance  as  a 

fun^t  on  of  w  .  The  sampling  plans  and  the  decision  rale  guarantee 

* 

that  this  probability  is  less  than  or  equal  to  1  -  P  for  all  W  £  w  . 

Suppose  there  exists  a  function  B(t,w)  such  that  F(t;w)  ^  B(t;W) 
for  all  t  and  all  w  ,  tnen  we  have 

(3.12)  P( Acceptance)  ^  L(B(t;w))  . 

In  order  t  insure  that  the  producer's  risk  will  not  exceed  a 
number  Q  (o  <  a  <  1),  whenever  w  ^  w,  ,  one  needs  to  satisfy  the 
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Inequality 


(J.l?)  L(B(t;0) )  £  L(B(t;W1))  £  1  -  Q 


which  will  be  satisfied  if  B(tjfl)  la  decreasing  in  0  .  One  may 

use  the  inequality  (3.13)  to  determine  c  such  that  for  fixed  ©^  , 

©^  and  P  ,  the  sampling  plan  and  the  acceptance  rule  will  insure 

that  the  producer's  risk  will  not  exceed  a  .  Alternatively  for 

* 

fixed  n  ,  t  ,  ©  ,  c  and  P  ,  one  may  want  to  know  the  minimum 
value  of  y  such  that  for  all  ©  >  ©^  ,  the  probability  of 

acceptance  £  1  -  a  . 

Ncv  we  describe  the  evaluation  of  the  X  function  for  the  IFR 
and  IFR  distributions. 

IFR  Distributions 

Case  (  i  )  ©  «=  Li 

r 


1  - 


e-t/<v 


l/r 


(MM 


B(t;nr) 


1 


t  £  H 


l/r 

r 


t  >  m. 


l/r 

r 


Thus,  for 


1  r 
p.  >  p  >  t 
r  r  c- 


ve  have 


(5.15)  L(F(t;ur))  £  L(B(t;pr))  £  L(B(t;^)) 


which  gives  a  lower  bound  on  the  OC  function.  This  bound  is  graphed 
beltjw . 
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and  the  1  ver  and  upper  bounds  can  be  obtained  similarly  by  using 
b(t;£^)  and  B(t;£  )  In  the  appropriate  intervals. 

Finally,  it  should  be  noted  that  sampling  plans  which  pr  )tect 
the  producer  can  be  constructed  for  test  time  less  than  the  g  al 
us .  ng  the  bound  B(t,y)  in  the  IFH  use. 

It  is  clear  fr  >m  the  b  unds  or.  UC  functions  f  r  trie  IFR  -ase 
that  we  have  to  test  beyonu  the  mean  quantile)  life  goal  t  protect 
tiie  c  nsumer  whereas  we  can  pr  tect  the  producer  by  testing  for  a 
time  less  than  trie  mean  \  quant,  le )  life  goal. 

DFR  Distr ibut.ons 

Case  i  )  *  =  u 

r 

In  this  case  there  is  no  nontrivial  lower  bound  for  the  uC 
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function  L  sir.ee  there  is  no  nontrivial  upper  bound  on  F  .  H.vever, 
there  existb  a  nontrivial  upper  bound  on  the  L  function  both  for 
p^  <  t  and  p.  ^  t  .  This  lower  bound  is  obtained  by  using  (3. 10) 

in  place  of  V  in  the  L  function. 

Case  ii )  *  =  T 

q 

Using  the  pper  l  jund 

1  ’  1  ^  Cq 

0.17)  B(t;C  )  =  |  t/C 

q  Ll-(l-q)  «  , 

one  can  construct  a  lever  bound  on  L  for  all  t  . 

Description  of  the  Tables  and  Seine  Illustrative  Examples 

Table  I.  This  table  gives  the  sample  sizes  necessary  to  estaolisn 

that  the  true  unknown  mean  p  of  an  IFR  distribution  is  at  least 

0  * 

p  at  a  nominal  confidence  level  P  .  Ihe  actual  level  attained 

is  also  tabulated. 

Illustration  1.  Suppose  an  experimenter  vants  to  choose  the  sample 

*  ,  0 

size  corresponding  to  c  =  2  ,  P  *  .  ,  >  •  t/p  •  l.u  ,  the.,  tius 

required  n  from  Table  1  Is  erual  to  5-  For  this  sampling  plan  tne 
actual  probability  level  Is  equal  to  .7523* 

Table  II.  This  table  gives  the  sample  sizes  necessary  to  establish 

th 

that  the  true  unknevn  q  quantltle  of  an  EFR  distribution  is  at 

u  * 

least  T  at  a  nominal  confidence  level  P  .  The  actual  level 

q 

attained  is  also  tabulated. 

Illustration  2.  Assuming  an  IFR  distribution  and  given  c  »  2  , 
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P*  -  .75  ,  X  -  t/C°  -  1.6  ,  q  -  .2  ,  we  find  from  Ifeble  II  that  the 
required  minimum  Maple  size  is  13.  For  this  sampling  plan  the  actual 
probability  attained  if  the  decision  rule  accepts  is  equal  to  .7961. 

Application  to  Drug  Screening 

Suppose  a  drug  which  is  known  to  be  non-toxic  is  administered 
to  patients  suffering  from  a  specified  cancer  to  determine  its  effective¬ 
ness.  The  drug  will  be  considered  effective  if,  say  20$  of  the  patients 
respond  to  treatment  after  T  weeks  (T  -  6,  say).  The  treatment  will 
continue  t  weeks  (t  »  say)  unless  more  than  c  patients  respond 
before  this  time --in  which  case  the  drug  will  be  declared  effective  at 
the  end  of  the  t  ^  weeks.  Let  F(t)  denote  the  probability  of  a 
response  by  time  t  and  F(T)  *  p  .  Since  we  do  not  want  to  discard 
a  potentially  useful  drug,  we  set 

c 

Plrejecting  drug  |p£  .20]  -  ^  (“)[F(t)]J[l  -  F(t)]n"J  £a  . 

j-0 

Nov  F(T)  £  .20  if  and  only  if  C  2o  ^  T  where  C  20  is  20th 
percentile.  If  we  assume  an  IFR  distribution  time;  i.e.,  j-  fffit) 
Increasing  In  t  ,  then  this  problem  is  identical  with  case  (ii) 
considered  under  IFR  sampling  plans. 

Illustration.  Ve  want  to  determine  the  number  of  patients  to  put  on 
trial  to  establish  with  90  per  cent  confidence  that,  F(T)  £  .20;  l.e., 

t 

the  response  rate  after  T  weeks  Is  at  least  20  per  cent.  If  ^  -  2  , 
c  «  3  ,  Q  -  .20  and  P#  »  .90  ,  we  find  from  Table  II  that  n  •  17 


and  the  true  confidence  level  Is  .91. 


Table  III.  This  table  gives  the  sample  sues  ..c  essnry  to  establish  that 


the  tru*»  nkr  -unntile  of  ft  .'FR  distribution  is  at  least  7 

q 

*  o 

at  a  nominal  confidence  J.evel  P  .  Here  X  «  t/q^  £  1  .  Also,  the 
actual  confidence  attained  is  tabulated. 


Approximations  for  Jample  Sizes 

The  sample  sizes  necessary  to  establish  a  quantile  for  the  IPR  and 
DFR  distributions  can  be  approximated  as  follcvs.  If  1  -  (1  -  q)^  is 
small,  then 


(3.10)  n  2  |  y-C-1-^ - t  ]  ♦  1 

'  L  1  -  (1  -  qr  J 
* 

where  yc+1  p*  is  the  P  percentage  point  at  a  standardized  gamma 

* 

variable  with  shape  parameter  r  «  c  ♦  1  or  one-half  times  the  P  per- 

2 

centage  point  of  a  x  with  (2c  ♦  2)  degrees  of  freedom.  This  approxi¬ 
mation  was  discussed  by  Gupta  (1962a)  where  tables  of  yc+p  p*  are  als° 
given. 

If  p  =  1  -  q  1  s  close  to  .5,  then  another  approximation  for  n 
is 


(3.19) 


“-  L 


(2c 


♦  1  ♦ 


Pip' 


)(l  ♦  J  1  - 


~2[ r 


(2c 

~P) 


1)2/(2c 


+  Pip* 


)2i 


J 


+  1 


* 

where  zp*  is  the  P  percentage  point  of  the  N  (0  ,  1)  random  varia- 

« 

ble.  For  example,  if  q=.l,X-l,c«l,P  ■•75,  then  *  = 

2.o92b  and  from  (3.18)  ,  n  -  27  .  The  exact  answer  from  Table  II  is  27. 
As  an  example  of  the  second  approximation,  let  X  «  1  ,  p  ■  1  -  q  -  .3  , 
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1 

i 


l 


* 

P  ■  .75  ,  c  «  2  ,  then  the  approximate  answer  from  is  1  which 

coincides  with  the  exact  value. 
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